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A cosmological gravitational wave background resulting from space-time quantum perturbations 
at energy scales of ~ 10 GeV is expected as a consequence of the general relativity theory in the 
context of the standard cosmological model. Initial conditions are determined during the inflationary 
(de Sitter) era, atz> 10 25 . A semi-analytic method was developed to evolve the system up to the 
present with no need of simplifying approximations as the thin-horizon (super-adiabatic) or the 
\Q , instantaneous transitions between the successive phases of domain of the different cosmic fluids. 

f"^ ■ The accuracy of such assumptions, broadly employed in the literature, is put in check. Since the 

physical nature of the fluid (known as dark energy) leading to the accelerated expansion observed in 
£S) , the recent Universe is still uncertain, four categories of models were analyzed: cosmological constant, 

X-fluid (phantom or not), generalized Chaplygin gas and (a parametric form of) quintessence. The 
results are conclusive with respect to the insensitivity of gravitational waves to dark energy, due to 
the recentness of its phase of domain [z ~ 1). The empirical counterparts of the gravitational wave 
, forecasts are still nonexistent for the noise levels and operational frequencies of the experiments 

^ — ■ already built are inadequate to detect those relics. Perspectives are more promising for space 

detectors (planned to be sensitive to amplitudes of ~ 10 -23 at 10 -3 — 1Hz). The cosmic microwave 
background is also discussed as an alternative of indirect detection and the energy density scale of 
inflation is constrained to be smaller than 10~ 10 in the analysis here presented. 

in 

00 

O 1 I. INTRODUCTION 

*sO ' 

The progresses accomplished in both experimental and observational fields during the last century permitted to 
2^ ■ establish a cosmological model in reasonable agreement with the reality (or at least with the glimpses provided by 
Oh| the experiments thereof). One of the most outstanding predictions emerged from that scenario is the existence 
of a cosmological background of gravitational waves. This forecast, almost as old as the general relativity theory 
itself (in 1916 Einstein published a work on linearized weak waves emitted by bodies with null self-gravitation and 
propagating through a flat space-time 0), was indirectly verified only in the seventies (ever since the binary system 
PSR 1913+16 has been monitored,showing an orbital deceleration rate compatible with a kinetic energy dissipation 
due to gravitational radiation 0, Q) an d up to the present no direct observational effort succeeded (a somewhat 
embarrassing fact in the so-called era of precision cosmology) . Nevertheless it would be hard to believe that this very 
consequence of the standard cosmological model (which explains the thermal history of the Universe, the genesis of 
large scale structures and the primordial nucleosynthesis) is incorrect. 

There is an immediate analogy between gravitational and electromagnetic waves. More important are the differences 
though. Only large massive objects in movement or coherent space-time vibrations emit gravitational waves, while 
electromagnetic radiation arises from incoherent superpositions of individual contributions of atoms, electrons and 
charged particles; near extreme gravitational fields (e.g., in the black hole surroundings) electromagnetic waves tend 
to darken, while their gravitational analog tends to be better emitted; unlike photons, that interact easily with 
matter, gravitons (the quanta of gravitational waves) interact very weakly, being able to cross regions as dense as 
the nucleus of a supernova or the primordial plasma just 10 -44 seconds after the Big Bang. These differences imply 
that if (not to say when) gravitational waves are detected and studied a completely new viewpoint (unaccessible to 
electromagnetic-only experiments) will be achieved and the consequences of that are unpredictable. 

The present work concerns the generation and evolution of primordial gravitational waves. Quantum perturbations 
of the space-time during a de Sitter inflationary phase are computed and taken as initial conditions for the subsequent 
evolution, governed by an excited oscillator-like equation. The evolution of those primordial perturbations up to the 
recent Universe is followed through a semi-analytic method. Since the gravitational wave background is stochastic 
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(isotropic and stationary) all the meaningful information is encapsulated in the frequency spectrum which is analyzed 
in sectionlVI Other authors |j,|5|,ISll3 accomplished similar calculations using analytic approaches which are compared 
with the results reported in section FVCl 

Perspectives of direct detection by both already built and planned experiments are treated in section IVI Al A 
promising indirect possibility also discussed in section IVIBI is the search for gravitational wave signatures on the 
angular spectrum of the cosmic microwave background (hereafter CMB), which corresponds to photons that decoupled 
~ 300000 years after the Big Bang, being sensitive to perturbations present in the cosmic plasma since that time. The 
analysis of possible observational counterparts characterizes the present work as a theoretical study able to establish 
forecasts for future experiments. 

About the recent Universe it is remarkable that the physical nature of the fluid (dark energy) representing ~ 70% of 
the present total energy density and leading to the accelerated expansion inferred by supernovae observations remains 
uncertain. A plethora of possible dark energy models have been proposed but the observational constraints are still 
not enough to determine its equation of state. Since dark energy seems to interact only via gravitation, a reasonable 
expectation could be to find signatures of its equation of state on the gravitational wave spectrum. The analysis for 
the most important models is done in section TVPI 



II. FORMAL DEVELOPMENTS FOR PRIMORDIAL TENSOR PERTURBATIONS 



In the linear regime, tensor (t) perturbations of a given flat, homogeneous and isotropic background (B) metric are 
represented by a transverse-traceless tensor hij to be added to the spatial part of the unperturbed metric tensor Q : 







(1) 



- where r, a(r) and 5tj are the conformal time, the scale factor and the Kronecker delta, respectively. In order to 
compute the evolution of these perturbations the Einstein-Hilbert action is expanded up to the second order : 

Som = ?nr* I d 4 xV^ l-g^hti.vhijj , (2) 



>0(2) = J " x Vtf v ' yy ' "/;.,."..;.!• 

neglecting the anisotropic stress tensor, which would act as a source term. The variational principle applied to 
implies that [llTol| 

h'l } + 2-h'ij - V 2 /i y = (3) 

(derivatives with respect to r being denoted by primes). One can also write J2J in terms of the cosmic time t = J a(r)d,T 
(with derivatives represented by dots) [Tl| : 

K ~ ~K - 2-hy - 4rV 2 h y - = 0. (4) 
a a a 2 

The conformal perturbation is related to its cosmic analogous by hy = a 2 (r)/iy and its Fourier expansion may be 
expressed as 

h£>(x,t) = VlG^G J h^(k,t) e W(k) expj-zk-x}, (5) 

where (A) accounts for the two polarizations, (+) and (x). For a wave propagating in the z direction, the polarization 
basis components are: 



= e x ® e y + e y ® e x , = e x ® e x + e y ® e y . (6) 

is syr 
.12] 



This basis i s sy mmetric, transverse, traceless and normalized - these properties are, respectively, expressed by the 
constraints 



e^(k)=e^(k), fe i ej ) (k)=0, e^(k) = and 



e W (k)e (V) «* (k) = 2(5aa/j e (A) (k) = e W* ( _ k) _ (?) 
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Using J5J), one may Fourier-transform J2J into 



h--h+ Po-2- h = 0, 



or into its conformal analogous 

h" + 2— ti + k 2 h = 0, h=h w (k,r) 



fc — — /x = 



or yet 



fi — h a. 



(8) 

(9) 
(10) 



Equation (|1(J|) above represents an oscillator excited by an effective potential a" /a and is a crucial point of the first 
simplifying assumption to be discussed in section TlV Al 

Considering once more the Fourier-expanded perturbation JSJ), the second order action J2J) is rewritten in a form 
which may be promptly quantized 01 : 



<? 0(2) = j2 1 — \h (X) 'hW*' - k 2 h^h^* 



drdk. 



(11) 



The conjugate momentum is = a 2 h^* , the corresponding operators are and Tr^. and the commutation 



relations to be satisfied are 



h W h (X "> 



"k ' ^k' 



iS XX 'd(k~k'). 



(12) 



Since h£ is hermitian its decomposition in terms of the creation and annihilation operators - considering the com- 
mutation relations |Q - is 

hi X \r) = h[ai X ^+h*[a^lr\ (13) 
where the evolution of h is governed by © and the standard commutation relations (for boson particles) are preserved: 

[«fe> V] = [«fe> a fe'] =0; l a k' a t'] = S kk'- (14) 

The power spectrum A 2 (fc, r) is easily obtained from the above formalism. Considering equations (|13|) and JSJ under 
the constraints imposed by the commutation relations (|14fl and the normalization properties J7J of the transverse- 
traceless tensor basis - and, just for sake of clarity, expliciting all the physical constants - the two-point correlation 
function is given by 



IMt,x)^(t,x')|0) = <Kx-x')£ 



V WnGh 
~2^f c 3 



2|/i (A) | 2 dk 



f&^Gh k 3 2 

<5(X-X') J -^- — \ h \ 2 dkLk, 



c 3 2tt 2 1 

where V is the volume of the system (normalized to the unity); the power spectrum is: 

2 d(0|^(r,x)^(r,x)|0) 
At(fc ' T) = dfofc 

and therefore, from l|15|l and (|16|l and recovering the natural system of units, it is straightforward to notice that 

A 2 (k lT )=64nG^\h\ 2 . 



(15) 



(16) 



(17) 



III. AN EFFECTIVE EQUATION OF STATE AND THE SCALE FACTOR EVOLUTION 



The standard (Big Bang + inflation) cosmological model states that (1) m the very beginning the Uni- 

verse experiences a de Sitter-like inflationary phase, when the perturbations are originated; soon afterwards, (2) the 
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inflationary scalar field decays quickly and the Universe becomes radiation-dominated (p r = l/3p r ); as the Universe 
expands, the radiation density falls as a~ 4 and at z ~ 1000 (3) the era of cold dark matter (CDM, p m — 0) starts; per- 
sisting until z ~ 1, when (4) the expansion becomes accelerated and the dominant fluid, dark energy, is characterized 
by a negative pressure. 

The first attempts to explain the dark energy tried to associate that negative pressure with the vacuum energy 
(A). The concordance model in that context would be a Big Bang + inflation + ACDM model. However, the 
concordance model faces great difficulties (the cosmological constant problem) and several alternative models have been 
proposed (see |l5l Il6| for review, and references therein). Here, four models were chosen for a comparative analysis: 
cosmological constant - the paradigmatic model, with an equation of state p = — p; X-fluid - a generalization over the 
A equation of state, p = cop, lo < 0, which is known as phantom fluid when u> < — 1 a first order parametrization 
(ui(z) = ujq + uj\z) of quintessence, which is a class of models where a scalar field is invoked to derive an evolving 
equation of state; and generalized Chaplygin gas, with an equation of state p = —A/p a = — Apo(po / p) a , A, a £ [0, 1], 
that behaves as matter for high z and as A in recent eras. 

Given the equation of state, the energy density p(x) and pressure p(x) evolutions are easily computed: 

P{x) = PioMx), P(x) = PMx), (18) 

i=r 1 m,de i—r^m : de 

where x = 1 + z = 1/a and the sum contains three terms: radiation (r), matter (m) and dark energy (de), this last 
one corresponding to one of the models mentioned above. The parametrization in terms of the functions <&i(x) and 
^i(x) is convenient to simplify the notation. With this parametrization, the Friedmann and Raychaudhuri equations 
assume, respectively, the form 

i i 

where H and fl io are the present values of the Hubble (H = a/a) and density (Qj = 3H 2 pi/8irG) parameters. The 
functions <&i(x) and ^i(x) may be written as: 

V r (x)=x 4 , $ r (x)=.T 4 /3, 

V m (x)=x 3 : <f> m (x) = Q, 

*a(s) = 1, M*0 = -1, 

9 x (x) = x 3 ^ +1 \ $ x (x) = m 3 ^ 1 ', ( 2 °) 

* gCg (x) = [A + (1 - A)x 3 ^] ^ , <$> gCg {x) = -A [* gCg (x)]- a , 

® q {x) = exp [3uji(x - 1)] x H^o-^i) j = [ Wq + CJl ( x - l)] y q ( x ) t 

where the sub-indexes refer to each of the state equations: radiation (r), matter (to), cosmological constant (A), 
X-fluid (X), Chaplygin gas (gCg) and quintessence (q). 

Therefore, the dynamics of the Universe may be computed with no need of truncating the equation of state in 
particular redshifts between a certain phase and the next one. The effective equation of state is defined as p = ui e f(x)p, 
where 

The evolution of u> e {(x) is plotted in figure ^ The present density parameter of each component is f2 ro = 5 x 10 -4 ; 
^m = 0.28; and £lde — 1 — ^m — ^Voi except for the Chaplygin gas, which acts as a unified dark matter-energy 
fluid: Sl g c ga = 1 — £lr - The dark energy parameters considered are: X-fluid: u> — —0.8; phantom: u> = —1.8; 
quintessence: ujq = — l,u>x = —0.2; and Chaplygin gas: A = 0.8, a = 1. These values are in agreement with the 
observational constraints presently available (e.g. see 0] and references therein). Figure Q also shows the step-like 
effective equation of state (hereafter step-eos, dashed line) which relates the second simplifying assumption to be 
discussed below. 

IV. SEMI-ANALYTIC APPROACH TO THE PROBLEM 
A. The thin-horizon and step-eos assumptions 

Equation Q has asymptotic solutions which are given by: 

r -i -1 

h(k, t) = exp(-ifer) a(r)v2k , k 2 > \a"/a\ r > r k (22) 
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FIG. 1: Effective equation of state in different cosmologies, all characterized by three fluids: radiation, matter and dark energy. 
The upper panel clearly shows the transition between the radiation- and matter-dominated phases and sketches the behavior due 
to the dark energy models investigated. The dashed line represents a model for which the phase transitions are instantaneous 
and the dark energy fluid is the cosmological constant. The lower panel shows in detail the effective equation of state at the low 
redshift zone of the upper panel for: cosmological constant (A), X-fluid (X), phantom (Xf), quintessence (q) and Chaplygin 
gas (gCg). 



h{k, T) = A k , k 2 < \a"/a\ r < r k (23) 

where A k is a constant and r k is the transition time when the k mode crosses the horizon {k 2 = \a"/a\). The thin- 
horizon (or super-adiabatic) approximation ignores the intermediate cases, k 2 ~ \a" /a\, joining the solutions (1221) and 



(12. 'it at r = rfe, by means of the condition A k = a(rfc)v2fc 



-l 



exp(-ikr k 



The calculation of cosmic gravitational waves using the thin-horizon approach is very simple - it is enough to 
compute the scale factor integrating l|19|l - and it becomes still simpler if the step-eos is assumed. These first 
calculations were accomplished by Grishchuk in the early seventies (see [TsL Il9j | for review) considering the radiation 
and matter phases only and have been recently extended to include a subsequent A phase 0- 

A formal treatment removing the thin-horizon assumption and taking an arbitrary number of successive step-eos 
phases (all with power law scale factor) was developed by Maia [2(j, who obtained the formal solution for 



ep-e 



h(k,r) ^h(k 1 T l )T(m+ 1/2) [kT/2] 1/2 - m J m _ 1/2 (fcr), m =^—^- (24) 

where is an initial time and fcr^ <C 1 and h'(k,Ti) = were assumed; T(m + 1/2) and J m -i/2(kT) are the Gamma 
and first type Bessel functions. Equation l|24|) is valid for equations of state p — up and does not include more 
exotic cases (as the gCg) . This formalism has been recently used |(| 0, EI to obtain the gravitational wave spectrum 
considering the ACDM phase. However, the step-eos is still maintained in these recent works. 
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B. Beyond the simplifying assumptions 



Using the fairly known relation x = l + z— ^,ao — 1, equation (jSJ) may be rewritten (performing a first change 
of variables from t to a, then substituting (|19f) into the resulting equation and finally applying a second change, from 
a to x) in the form: 



dx 2 



3x 2 A\dh 
~2~~B J dx~ 



2 

:r 2 



3x- 



A 



V _ 

, B ~x 2 "' B 

i—m,r,ee 

where k is the (non-dimensional) wave number: k = fc/i?Q. A convenient new variable 77 may be introduced taking 



i—m,r,ee 



h = , 



(25) 



hx 2 = 77 exp 



as a definition. The resulting equation is similar to (|10|l 

d 2 V 



3Ax 2 /4B) dz 



dx 2 

A = 2/x 2 + fx 4 + f-x 2 
dx 



(k 2 /B + A) 77 = 0, where 
B = B, f = A/B. 



(26) 



(27) 



Therefore, solving equation l|25l) is equivalent to integrate l|27[) using l|2t)[l to recover the solution for h. Formally, 
(|27H and (|10|l are similar to each other and at a first glimpse this new formulation offers no advantage over that one. 
However, a second look at (|27|) shows a different scenario. 
For an effective equation of state p = u) fP, one has 



k 2 /B + A= k 2 x 3n + 3 + (n + 2)(n + 3) X - 2 , 

If n does not vary with x, then the solution of (|27|1 is analogous to (|24|l : 

1 



71 = — W e £ — 2. 



(28) 



77(1) = 




+ c 2 r(l + /)J; 



' 5 + 3™ 

2kx~ ~ 
5 + 3n 



(29) 



. V23 + 20n + An 2 

1 

5 + 3n 



The assumption u) e t = constant is valid for sufficiently small intervals of x and a feedback process may be used to 
consider its variations. Establishing a set of solutions such that each element ?7j(x) corresponds to an interval [xj, Xi+i], 
the solutions are all given by (|29|l but with varying values of (ei,C2)j. A set of constants (c\,C2)i is associated to 
the solution set. Taking the initial conditions from a given inflationary model, the constants (01,02)0 of 770 (x) are 
obtained. This first solution is valid since Xo = Xk, when the k mode crosses-out the horizon. A recurrence relation 
schematically expressed as 



»7i(Xt) = T)i-i(Xi) 

n'iixi) = V'i-i(xi) 
Xi — Xi-i + Ax 



> (ci,c 2 )i 

gives the remaining constants and synthesizes the semi-analytic method here employed. 



(30) 



1. Initial conditions 

The necessary initial conditions are defined in the inflationary era (see e.g. |22T | for review), when the super-adiabatic 
approach is applicable. During the slow-roll phase one has H w and in the exact de Sitter case, H — 0, so that the 
asymptotic regimes k 2 <C \a" /a\ and k 2 ^> |a"/ a l are equivalent to k <§; aH and k 3> aH, respectively. 
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Assuming that the perturbations start all inside the horizon - with physical wavelength smaller than the Hubble 
radius (k > aH) - it is possible to take the solution IL'l'li for r > t&. Ignoring the oscillatory part, the initial condition 
to be taken is 

hk{Tk) = ^tm = ^T^ k = a(r k )H, (31) 

where a(r) and, consequently, H depend on the inflationary model adopted. In the slow-roll regime H is proportional 
to the potential of the inflationary scalar field and for the de Sitter model this potential is exactly constant so that 



H 2 = ^V(^ = ^Mm% (32) 



where M is the energy density scale of inflation. The energy scale E v is the only free parameter of this model for 
which V(<p) = and E v = M 1 / 4 m P i. 

Equations (|31|l and i|32|) lead immediately to the result 



hk(r k ) = \ —— m Pl k-*' 2 & W{x k ) = J-r- mpi k- 3 ' 2 x^ 2 . (33) 



k 

h' k (r k ) = * K(x k ) = -yj^- mn ^ m ~~ 3 



x- 3 . (34) 



To conclude the initial conditions set-up an expression for Xk is needed, for each mode k crosses the horizon in a 
different redshift. That is done using the condition k = a{x^)H / Hq, so that 



Xk 



y-g- mpi {H k)-\ (35) 



V. GRAVITATIONAL WAVE SPECTRUM: RESULTS AND DISCUSSION 
A. The power spectrum redefined 

The power spectrum (fT7|) is calculated using the semi-analytic method described above: 

A 2 (k 7 x) = ^f(H k) 3 xt\x 2 K(x)\ 2 . (36) 

Knowing that the largest observable wavelength corresponds to the size of the current horizon, 1/Ho, the lower limit 
of physically interesting frequencies is obtained: k = 1 or v = Ho « 2 x 10~ 18 Hz (this was the first k leaving the 
horizon and did not return until 2 = 0). The upper limit, in turn, corresponds to the last mode leaving the horizon 
at the end of inflation. Considering that the radiation era starts immediately after inflation (instantaneous and 100% 
efficient reheating processes) and that p r = /O ro a; 4 , the end of inflation is: x en d = (87rM/3)mp;/(i?o^ro)- Using l|35|l . 
the maximum wave number is immediately obtained. 

Another important quantity is the amplitude, correlated to A|(k, x) through: 



h c (k,*) = ^ifel (37) 

and of great interest when the detectability of primordial gravitational waves is computed. The amplitude is ap- 
proximately equal to the variation observed in the distance L between two proof masses submitted to a gravitational 
radiation flux: h c (k,x) w AL/L. 

The third quantity of interest is the spectral energy density, defined as 

3H 2 d(0\ Pgw \0) 

n gw (k,r) = — ^ , Pgw = -T . (38) 

To calculate Cl gw (k,T) the action (|llfl is considered to obtain the stress-energy tensor and the same procedure that 
led to tHHl can be applied to obtain, 



(0|p 9W |0) = 



fc3 (/^.) 2 + fcXdfc 







2tt 2 a? 

o ( h T) - ^ G fc3 W + *X _ 1 k 2 Al(k,r) 
9w[ ' ' ~ ZH 2 2^ a 2 ~12 a 2 tf 2 ' 
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where in the last term of (|39|l it was considered that \h'(k, r)| 2 = k 2 \h(k, r)| 2 . As a function of k and x, the spectral 
energy density is 



Qgw (k, &) 



x 2 (# k) 2 A 2 (k,x) 



12 



H 2 (x) 



(40) 



The spectra thus defined are all interesting. They contain information on the amplitude of the tensor perturbations 
, their detectability ljo7|) and the energy density of each mode lj4TJ|) . 

B. During the inflationary era 

According to the initial conditions (|35|l and l|o6|l . the amplitude spectrum at x — x cnd is flat, as shown in figure [21 
for four energy scales: M = 1CT 24 (E in{ = lO- 6 m P i); M = 1CT 20 (E in{ = 10- 5 m P i); M = 1CT 16 (E inf = 10- 4 m P; ); 
M = 10~ 14 (E inf = 10~ 3 - 5 m P; ). The spectra were calculated at x cnd (M) » 10 25 - 10 28 . 

Defining <5 h = Alogh c /logh c and 5m = AlogM/logM, one notices that 8^ « 3/4<5j\/. For x < x cnd , the spectrum 
evolves such that the amplitude of modes inside the horizon is always decreasing (sections IV CI and IVDfl , while the 
external ones remain almost constant. This fact allows, in principle, to establish an upper limit for M. Since the 
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FIG. 2: Amplitude spectrum of relic gravitational waves taken immediately after inflation for four energy scales: M = 10 
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metric perturbations leave their imprint on the CMB anisotropy spectrum, the maximum amplitude should not be 
larger than AT/T « 10~ 5 at x ~ 10 3 — 10 4 , when photons decouple. This constraint implies that M < 1CP 14 (this 
conservative estimate is discussed in section lYl Bfl . 

A remarkable characteristic to be noticed in figure [5] and the following ones is the oscillatory shape of the spectra. 
Far from being a purely numeric effect, those oscillations result form the very mathematical solution of ©: there are 
oscillations not only in time, but also in the Fourier space. 



C. During the non-accelerated expansion 

The result of the semi-analytic procedure is shown in figure |3 for different redshifts inside the radiation era. The 
modes of smaller wavelength cross-back the horizon at larger redshifts so that the modifications in the power spectrum 
start at the higher frequencies and the information about the initial conditions are better preserved for low frequencies. 
In the area corresponding to the modes that already crossed the horizon, the slope is a k, —1.2. 

The transition from radiation- to matter-dominated era begins at about x « 10 and finishes at x « 10 2 (figure^). 
It is important to highlight that the shape of the spectrum carries information mainly about: (1) the initial conditions 
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Log v Hz 




Log k 

FIG. 3: Tensor power spectrum (chirp amplitude), for M = 10~ 14 , at different redshifts inside the radiation-dominated phase. 
The slope of the spectrum is a ~ —1.2, which is comparable to previous estimates 

determined in in the primordial Universe and (2) the evolution of the effective equation of state - or equivalently of 
the scale factor. While out of the horizon each k is maintained (almost) unaffected preserving its initial value; when 
returning to the horizon, it begins to evolve; however, when l/H is much larger than k, l|29l) tends to the asymptotic 
regime (|23|l . when all of the modes fall proportionally to a^ 1 . For this reason, one observes that the high frequency 
region in figure [21 maintains the same slope (only the amplitude varies) during the subsequent phases. 

Figure 0] shows in detail the low frequency region of the power spectrum calculated at x — 5, far inside the matter- 
dominated era, when almost the whole spectrum has been modified from its initial condition. The spectrum has an 
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FIG. 4: Low frequency region of the amplitude spectrum of primordial gravitational waves, calculated at x = 5, for M = 10~ 14 . 
The dashed line corresponds to the result expected from the step-eos assumption. 

almost null slope in the region 10 1 < k < I0 24 while in I0 24 < k < 10 4 its inclination is modified because of the 
variation of u> e {. That feature is completely different of what is obtained from an step-eos. Modelling w e f with a 
step-function (dashed line in figure^), other authors |EIEEi0| obtain a spectrum with a ps —1.5 for k < I0 3 [||. 
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For comparison a dashed line is included in figure corresponding to the slope obtained with the step-eos. 

A last comment is concerned with the region 1 < k < 10 of figure 0] The picture was taken in x — 5 and so the 
modes within that region are in the threshold of the transition from "outside" to "inside" the horizon (k ~ aH). Since 
the thin- horizon simplification is not considered in this work, that region shows a transitory slope (a rs —1.6), that 
disappears for x smaller than 5, as the modes enter the horizon. This characteristic is not perceptible in the previous 
illustration due to the scale. 



D. During the accelerated expansion 

Nearly all the modes are already inside the horizon when the accelerated expansion begins. For that reason, the 
shape of the spectrum remains constant and only the amplitude is altered in that phase. Figure El shows a narrow 
slice of the current spectrum, so that the effect of the different models of dark energy can be better observed. 
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FIG. 5: Present spectrum [x = 1) of gravitational waves for M = 1(T 14 . The dark energy models considered are: cosmological 
constant (solid), Chaplygin gas (blue), X- fluid (green), phantom (dashed), and quintessence (gray). The parameters used are 
the same as in figure Q 

The first characteristic to stress with regard to figure is the complete degeneracy between cosmological constant 
and quintessence models (u>q = — = —0.2), whose curves overlap completely. Moreover, the Chaplygin gas 
(A — 0.8 and a = 1) is also almost degenerate with those two curves. Observing figure it is easy to notice that 
these three cases have very similar behavior and this degeneracy is expected. The X-fluid (u> — —0.8), in turn, 
leads to an amplitude smaller than the cosmological constant while the phantom fluid (lj — —1.8) acquires a larger 
amplitude. The difference is greater for phantom than non-phantom fluid because the chosen parameters are such 
that \u>x t ~ w a| > l^x — wa|. Since the dark energy begins to act very recently x ~ 2, the signature of different 
models is very weak - to produce a difference of « 2 in log h c scale, it is necessary to vary uj of almost one unit (the 
more negative is the equation of state, the larger is the amplitude). 

Although just a particular choice of parameters is shown, other possible combinations have been tested. The results 
allow to state that: (1) for gCg l24j the parameter a does not have any visible effect, while A increases the 
amplitude as it varies from (corresponding to the minimum stipulated by the spectrum without dark energy) to 1 
(reaching the ACDM curve); (2) the X-fluid [2{| has only one parameter, ui, that plays the same role of A, but it 
assumes values more negatives than —1 and, in that case, the generated spectrum is above ACDM; (3) quintessence, 
in general, is completely degenerate with the X-fluid, unless that uj% ~ 1 (or greater) is assumed, but that does not 
have any physical sense for in this case the quintessence would affect the dynamics of the early Universe, being in 
contradiction with the most accepted cosmological models. 



11 



VI. SEARCH FOR EMPIRICAL COUNTERPARTS 



A. Direct attempts 



The term "direct detection" should be strictly used only for experiments whose principle is to measure the deforma- 
tion of the space-time submitted to tensor perturbations. The interfcromctric experiments lay in that class. However, 
it is common to include under that label detectors that measure the deformation of a proof mass in resonance and 
that is done in this work. The quantities A 2 , fl gw and h c characterize, in an equivalent way, the gravitational wave 
background and do not depend on the detector. However, the signal S(t) registered by an hypothetical detector will 
have a component s(t), representing the observable, and another, the noise r(t), which depends on the characteristics 
of the detector: S(t) — s(t) + r(t). The mean (square) contribution of these components are |26|: 



p r°° r°° 
{s 2 {t)) = - / u-Xdv, (r 2 (t)) = / h*du, (41) 

z Jo Jo 

where F is an efficiency factor representing the loss of sensitivity due to the fact that gravitational waves come from 
all directions, while the detector is only maximally sensitive in some preferential ones. That expression for (s 2 (t)} is 
valid for a stochastic background and the gravitational radiation may be considered detectable if 

h c > Jyh r . (42) 



1. Resonant masses 



Schematically, a resonant mass detector is a (cylindrical or spherical) massive solid mass whose mechanical oscil- 
lations are excited by gravitational radiation and converted to electric sign by a sensor. The variation of the mass 
length is a sum of all its vibration modes, but the sensor is filtered to receive only the fundamental frequency, being 
therefore a narrowband detector. 

In the case of a cylindrical bar, the efficiency factor is F — 8/15, for the preferential direction of incidence is any 
direction perpendicular to the bar (if the sensor is located at one of its ends). The sensibility of the last generation 
experiments is^h r fa KT 22 Hz~ 1/2 , at ~ 900Hz. Using g2Jl the minimum detectable amplitude of those experiments 
is: 

h c > 5,8 x 1CT 21 , for k « 4 x 10 20 . (43) 

A spherical detector has a larger mass for the same resonance frequency implying in a greater cross-section. Besides, 
it is sensitive to different directions and polarizations of the incident radiation. The sensibility of those detectors is 
planned to be j^j h r ~ 10~ 23 Hz -1 ^ 2 , for frequencies between 200Hz and 2kHz. Equation l|42|) here implies that 
resonant spheres are able to detect gravitational waves of up to 

h c > 2, 7 x 1CT 22 , for k w 1 x 10 20 . (44) 

Comparing bars 1|43|) and spheres (|44(1 . one notices an improvement of one order of magnitude, but that is still not 
enough to detect waves of cosmological origin. 



2. Interferometers 

Laser interferometry is used in large detectors to measure the displacement of free-falling masses submitted to 
gravitational radiation. The basic components of those detectors are two arms of length L connected at one end and 
three masses (one at the junction point, two at the free ends), constituting three pendulums. The preferential incident 
direction is perpendicular to the plan of the arms and F = 2/5 sin 2 9, where 9 is the angle formed by the arms [2f|. 

Unlike resonant masses, interferometers are broadband apparatus whose maximum wavelength is limited by the 
size of the arms. The main interferometric detectors are: LIGO j2(|, with a 4km arm length; VIRGO [30l l3l| . 3km- 
sized; GEO600 H2 and TAMA300 HH, with 600 and 300 meters, respectively; AIGO/ACIGA [13, planned to have 
4km; and finally, the most ambitious projects, BBO and LISA j^| are space interferometers where three arms 
(5 x 10 6 km in LISA and 5 x 10 5 km in BBO) are arranged in an equilateral triangle. 

The detectability estimate is performed using (|42|) once more. The results are shown in figure [S] In the high 
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FIG. 6: Gravitational wave spectrum for three inflationary energy scales (M = 10 1 , 10 1 , 10 1 , from bottom to top, thick 
lines) in comparison with sensitive estimates for the detectors. 

frequency band, resonant masses and ground-based interferometers are shown, all with insufficient sensibility to 
detect primordial gravitational waves. The sensitivity curves of resonant masses are nearly reduced to a point in that 
scale, for these are, by construction, narrowband detectors. The interferometric detectors are also identified in the 
figure. Only LISA and BBO would be able to detect a cosmological signal, but both are still in project [35ll36T| . 



B. CMB: an indirect attempt 

In simple terms, one can say that |37| until « 10 3 , the temperature of the Universe is high enough to ionize 
the hydrogen and, via Compton scattering, photons and electrons are coupled; in addition, electrons are also coupled 
to barions through electromagnetic interaction. The radiation pressure offers resistance to gravitational forces and 
acoustic oscillations arise in the plasma. In z*, the hydrogen recombines (and the photons last scatter); the Universe 
becomes then transparent to photons and compression and rarefaction areas of the plasma at that redshift represent 
hot and cold areas, respectively; besides, photons suffer gravitational redshift when leaving the potential wells of the 
last scattering surface (Sachs- Wolfe effect H3). The resulting fluctuations appear as (primary) anisotropics in the 
sky. 

The theoretical calculation of the CMB anisotropies is based on the linear theory of cosmological perturbations [l(| • 
The temperature anisotropy at position x and direction n, Ay(x, n), depends, in principle, on both the direction and 
the frequency, but the distortions in frequency are of second order. A Fourier expansion of At(x, n), results in modes 
that propagate independently of each other. Assuming axial symmetry around k, a Legendre expansion can also be 
done: 

A T (k,n) = ^(2Z + l)H) i A T! P / ( M ), /i = k-n/fc, (45) 

i 

where Pi(/J,) is the Legendre polynomial of order I, At, is the associated multipole moment and A*r(k, n) is the Fourier 
transform of Ay(x, n). A similar expression can be written for polarization anisotropies Ap(k,n). 

The Boltzmann equation, which describes the temporal evolution of the Stokes parameters of the radiation field, is 
constituted by a collisional (referring to the Thompson scattering) and a non-collisional term. In the case of tensor 
perturbations, the resulting equation is (3^ 



A'^ + ikfxA^ = -ti - k' (A$) - Y) 
A'f + ikn&$ = -«' (A^ + y) 



(46) 
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Y = 

where h is governed by equation (jSJ) and 



!a® + -A« + J-A® - * Ag> + -A« - J-Ag> 
10 To 35 T2 210 T4 5 p ° 35 Pa 210 p " 



TO 

K'dr (47) 



represents the optical depth (and k' = an e x e <TT)- The solutions of (|46J) may be written in the form [40T | 

/■TO 

T{k,r )= e tk ^ T - T ^S(k,T)dT. (48) 
Jo 

where J 7 is the (temperature or polarization) anisotropy, To is the present time and S is the corresponding source 
term: 

= -e- K h' + ICY S { p } = -ICY, (49) 
where JC = «/exp(— k) is the visibility function. Expanding T in terms of tensor multipoles one has 

Fl(k,ro) = £ y|^T t fc ( T - T oT 2 ji[k(r - T )}S(k,r)dT . (50) 
If the observation takes place at x = 0, the expansion of T in spherical harmonics is 

A. i 1 /" 

f = EE a lm Y lm (6, 0) = — — / dkF z , (51) 

j[=0 m— — l ^ ' 

where Y\ m is the angular part of the spherical harmonic eigenfunctions, such that: 

Pi(j*) = X^l52Y£ n (e,<t>)Y lm (e,<j>), [i = cos8 (52) 

TCI — — l 

^From lp)2"|) and (|FT]l it results 

fl ' m = " (27t) 3 / 2 y ^ k ^ ^/m(^)^)> ( a iniai'm') = Ci8w5 mm > ■ (53) 

The evolution of the Boltzmann equations does not depend on k, so it is possible to take , ?i — * h^Ti , where 
h m i is the initial condition of the perturbation. The correlation function is (ft.i n i(k)ft,* ni (k')) = A 2 (k)S(k — k') and so 
the Ci coefficients arc determined: 



47T 



(2tt) 3 / 2 



J A 2 |^| 2 (Y^Y^dk = ^J A 2 |^| 2 fc 2 dfc (54) 



In short, to calculate the induced tensor perturbations of the temperature (polarization) Ci spectrum it is necessary 
to compute the multipoles Ay (Ap') using (|50|) and the source term 5^ («Sp ); the result is applied directly on 154|l . 
using the initial power spectrum A 2 , given by equation (|17|) and shown in figure |21 All information on the evolution 
of gravitational waves after inflation is contained in the source terms. 

From 149|) it is noticed that the polarization source term is suppressed if the media is optically thin, so that after 
reionization the polarization due to gravitational waves is negligible and the polarization spectrum may be used to 
verify the initial conditions obtained in section IIVB II Although both tensor and scalar perturbations affect the 
CMB, it is possible to find a mode (B polarization mode) which is induced only by gravitational waves and this is a 
promising future perspective to constraint the initial conditions of the problem. 

The temperature source term in 1)490 is reduced to <S™ ~ —h' after recombination and, therefore, it contains 
information about the evolution of the power spectrum since then. Moreover, after the condition k <C aH is achieved 
the amplitude h falls as a -1 and only the modes k < fc* (where fc* is the wavelength of the order of the Hubble radius 
at recombination time) contribute to the source term . 
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The power spectrum of temperature anisotropies was calculated according to the prescription above, using the 
profiles of k and K. obtained from the Saha equilibrium equation and the envelope of the spectrum corresponding to the 
ACDM model with M — 1CP 14 . The integration limits (fc m i n = 1 and fc max = 1500 ~ fc», x m - m = 1 and x max = 1500) 
were chosen in order to reduce the computacional cost without losing relevant information. Only the multipolcs 
I < 10 were calculated and the result is a plane spectrum, with amplitude Cil(l + l)/27r = 41.1(/iK) 2 . The form of the 
spectrum is compatible with other theoretical predictions 0, , as well as with observational results > but the 
value found is ~ 25 times lower. This means that the inflationary energy scale M could be larger than the conservative 
limit adopted along this work. With M — 10 -10 , this difference is suppressed: Cilil + 1)/2tt = 1.01 x 10 3 (/xK) 2 . The 
current spectrum of gravitational waves can have in this case the amplitude given by the superior curve of figure 

The limit M — 10 _1 (E^ = 10~ 3 ' 5 mpi) is a direct consequence of the condition that the amplitude of tensor 
perturbations should not be larger than the observed CMB anisotropies and, at a first glimpse it may seem in 
contradiction with the result presented above. Under a more careful glance, however, one notices that this last result 
is obtained when considering the physical processes involved in the generation of anisotropies induced by gravitational 
waves. Additional physical ingredients (specifically the ones regarding to the recombination process) that determine 
the efficiency with which the perturbations imprint themselves on the CMB spectrum had to be considered. In order to 
avoid the additional uncertainties introduced by these new ingredients, the authors have chosen the most conservative 
limit, M — 10~ 14 , when performing the analyzes of section IVBI However, other possible values have been considered 
in the literature between M ~ 10" " 8 and M ~ 10~ n 7], which are consistent with the results here presented. 



VII. SUMMARY AND CONCLUSIONS 



The origin and evolution of the primordial gravitational wave background was computed for different cosmologies. 
Temporal evolution of these relics is governed by an oscillator-like equation with an exciting potential a" /a. The 
initial conditions are set-up during inflation and, for a de Sitter model, the initial spectrum is flat. As the modes 
cross back to the horizon, ever since the end of inflation until x ~ 10 4 , the spectrum assumes a slope (a rs —1-2) in 
the high frequency region (k > 10 4 ). That slope is maintained in the subsequent eras. The modes with k between 
~ 10 4 and ~ 10 2 (a = —0.75) enter the horizon during the radiation-matter transition and the step-eos assumption 
is shown to be unsupported in these redshifts. The final spectrum is flat in the low frequency band 1 < k < 100. 

The analysis here accomplished also shows that gravitational waves are not strongly affected by dark energy. 
Establishing constraints on dark energy using gravitational waves is not a promising task (even though one could 
arrive to the very opposite conclusion naively considering that both of them interact only via gravitation) . One must 
measure the gravitational wave spectrum with precision better than 5i O g 10 h c J$ 0,4% to constrain u> within 8 U < 1% 
(with the energy scale M independently fixed and an X-fluid equation of state assumed a priori). 

Once obtained the theoretical spectrum of cosmic gravitational waves, the natural question concerns the empiric 
counterparts of those forecasts. None among all the already built experiments is sensitive enough to detect primordial 
gravitational waves, because at their operational frequencies that signal is orders of magnitude below their noise levels. 
The perspectives are promising for the future space detectors though. 

CMB, in turn, supplies a concrete possibility of obtaining this information indirectly, but the influence of gravita- 
tional waves on the CMB must be understood in detail in order that one can be able to extract information from the 
angular spectra whose detection is expected to be greatly improved in the near future. The calculation here performed 
establishes a limit on the energy scale of inflation: M < 10 -10 . 

This work was financially supported by FAPESP and CNPq. 
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